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Abstract 

A common problem in the optimization of structures 
is the handling of uncertainties in the parameters. If 
the parameters appear in the constraints, the uncer- 
tainties can lead to an infinite number of constraints. 
Usually the constraints have to be approximated by 
finite expressions to generate a computable problem. 
Here, using the example of the topology optimization 
of a truss, a method is proposed to deal with such un- 
certainties by using robust optimization techniques, 
leading to an approach without the necessity of any 
approximation. With adequately chosen load cases, 
the final expression is equivalent to the multiple load 
case. Simple numerical examples of typical problems 
illustrate the application of the method. 

1 Introduction 

Since the first work by Michell [1904], a great num- 
ber of methods for the optimization of truss geome- 
tries were developed. A good overview of methods is 
given in Bendsoe et al. [1994]. A possible way to deal 
with uncertainties in parameters during the optimiza- 
tion process is robust optimization, a detailed survey 
of which can be found in ^eyer and Sendhoff [2007]. 
The idea behind the approach chosen in this paper for 
the robust topology optimization of a truss — follow- 
ing the ground structure method [Dorn et al., 1964] 
with the use of a matrix force method [Przemieniecki, 
1968] — is outlined in the following. 

A general optimization problem mmxz{x) subject 
to X G S is often given for various different param- 
eters, and principally represents the structure of the 
problem. For example, the sizing optimization of a 
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truss determines the overall structure. The meaning 
of the structure arises from the specification of actual 
data, for instance the topology of the bars, the defini- 
tion of the supports as well as the loads. In general, 
one has the optimization problems mirxx z{x,y) sub- 
ject to X E S(y) for different fixed parameters y de- 
fined by the actual data. 

If the parameters are not exactly known, one has to 
consider min^ z(x,y) subject to x G S(y) for y G Y. 
Here, Y represents the set of the uncertain parameters. 
Even without stating precisely the origin of this un- 
certainties (or scattering) in the parameters, the ques- 
tion arises on how to obtain a computable optimiza- 
tion problem. 

The reasons for uncertainties in the parameters can 
be e.g. measurement errors if the parameters are actu- 
ally obtained as or derived from measured variables. 
A parameter could also scatter across intervals which 
are restricted by limits provided by the user. Later it 
will be shown that it can make sense to choose an ar- 
tificial limit for the parameters. 

The importance of uncertainties arises simply due 
to their occurrence as e.g. measurement errors in real 
live applications. Measurement errors can basically 
be divided into two groups: systematic and random 
errors [Taylor, 1997]. Systematic errors are caused for 
instance by bad calibration of instruments or defective 
design. The distinctive feature of this kind of error is 
that it is always directed, and its impact is not nec- 
essarily known, which makes it difficult to treat in the 
error analysis (a miscalibration can cause an offset of a 
measured parameter, and this is not a scattering), and 
the best strategy is to identify and remove all sources 
for systematic errors. 

Random errors originate from a multitude of 
sources. Common sources are the thermal noise su- 
perimposed on a measurement due to a finite mea- 
surement device temperature, statistical fluctuations 
in counting experiments, or limited resolution of in- 
struments yielding upper and lower limits for a pa- 
rameter, to name a few examples. Random errors can 
usually not be avoided, but if their nature is known, 
they can be treated in a defined manner. 

If uncertainties follow a known statistical distribu- 
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tion, this knowledge can be integrated in procedures 
handling the parameters. For example in the case of 
the calculation of a mean value from a large number of 
independent measurements, a Gaussian distribution 
is assumed to justify the average as the most probable 
value. Other well-known and common distributions 
are the Poisson distribution in counting experiments, 
or the equipartition of the dice roll. 

If the distribution is completely unknown due to 
the lack of a justifiable model, boundaries of the in- 
terval in which the measured variable lies have to be 
established to generate a treatable problem. This is the 
equivalent to a user supplied limit for a parameter. 

In the case of an either unknown distribution or if 
only the limits are known, one could choose an arbi- 
trary nominal value yo G Y and solve the problem for 
it. This procedure can lead to results which are in- 
valid for any other value drawn from Y. More rea- 
sonable is the consideration of all parameter values 
y E Y - matching the worst-case-approach. This im- 
mediately leads to the robust optimization [c»en-ia. 
and Nemirovski, 19^ ] constituting the main body of 
the work presented here. 

Even in the case of a known stochastic distribution, 
a robust optimization can be carried out by the in- 
terval boundaries if required, although then the in- 
formation of the distribution is neglected, which can 
be considered by more reasonable procedures. The 
expectation value can be used to specifically choose 
a nominal value and solve the problem with respect 
to that value, but again accompanied by the loss of 
a large part of the available information. Due to the 
characteristics of the expectation value the approach 
considered here constitutes an average-case; however, 
the expectation value is just one special value, and 
therefore possibly leads to a solution valid for one 
point only. An example for the difficulty of this single- 
value approach is given by the problem of construct- 
ing a building at the coast under the assumption of 
a typical (ever-present) wind load (e.g. the expecta- 
tion value), but neglecting the possibility of instanta- 
neous absence of wind, causing the building to col- 
lapse if the latter case occurs. Another possibility 
is to consider only some very probable cases (larger 
than same given probability po) and thus to disre- 
gard the improbable cases. For very small po this is 
called reliability based design optimization (RBDO), 
while for larger po the denomination is robust design 
[Doltsinis and Kang, 2C ]. Formally, this leads to 
miuj z(x,y) s.t. x E ^p{y)>pQS{y). The robust opti- 
mization is again a direct implication. Nonetheless 
one should not neglect the sometimes serious con- 
sequences of unlikely cases, even if they seem to be 
statistically improbable in theory. Therefore, this ap- 
proach does not in general lead to safe-life, fail-safe or 




Figure 1: (a) Ground structure with fixed nodes on the 
left side and arrow representing the force (example 1). 
(b) Optimal solution. 

damage tolerance. 

Another possibility for the source of uncertainties 
are application errors by the user. With the assump- 
tion of a range or an order of magnitude of these un- 
certainties this case can be treated as above and leads 
again to robust optimization. 

First, we will look at the expected numerical be- 
havior (condition) of topology optimization of truss 
[Przemieniecki, 1968] in view of the topology to ob- 
tain artificial limits for the robust optimization. In (3)- 
(6) we will present the typical topology optimization, 
which will be extended to a multiple load case in (10)- 
(13) by a — for the robust optimization — unneces- 
sary approximation. We then start over with (14)-(17) 
as the new initial problem to reach its robust counter- 
part (18)-(21), which will be solved without approx- 
imation. The robust counterpart (18)-(21) is a semi- 
infinite optimization problem. By looking at the finite 
numbers of design variables, namely the cross sec- 
tions of the bars, we will be able to bypass the infinite 
numbers of constraints and obtain a linear program 
(32)-(35), which provides a worst-case or robust solu- 
tion. 

2 Condition of Topology Opti- 
mization of Truss 

The simple example of a 2-dimensional topology op- 
timization of a plain frame (cf. (3)-(6)) on the ground 

structure of the 3 bars Vj = (^4^^^^^^ with the nodes 

ej^) = {OJ)T e ]R2 and = (1,2)^ G IR^ for 

(1) 

i = 1,2,3, which are fixed in , / = 1,2,3, demon- 
strates the bad condition concerning a perturbation of 
the force direction in the case of a horizontal force act- 
ing on the node (1,2)^ with regard to the topology 
(Fig. 1(a)). 
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Figure 2: (a) Optimal solution of example 1 (c.f. Fig. 
1(a)) for the case with a perturbed force, (b) Robust 
optimal solution of example 1. 

For the case of an unperturbed horizontal force of 
quantity / G IR acting in positive x-direction, the opti- 
mal solution, illustrated in Fig. 1(b), yields the thick- 
ness of the bars as 

sW=(0 X o)', (1) 

where a- < < represents the stress limits 
(e.g. lower and upper yield point), which is deter- 
mined equally for all bars. For a perturbation in pos- 
itive y-direction, the force / is now constructed with 
the components fx > fy) the magnitude of the force 
\\{ fx fy )\ \ — f remains unchanged. As one opti- 
mal solution ifor — (7_ < C-^ the thickness of the bars 
results in 

§(2)^ / ^ M )^ (2) 

which represents another topology than the unper- 
turbed case (Fig. 2(a)). Since the optimization 
problem is a linear program, there are typically 
many solutions with the same objective function 
value. The stated solution s*^^) is an edge of the 
feasible region. Obviously, another edge is s^^^^ = 

^ Q ^ and every point between these 

solutions is a solution, too. 

Apparently, here the solution of the unperturbed 
case is not feasible for the perturbed case. The naming 
convention of the perturbed and unperturbed case is 
of philosophical nature. Therefore, a numerically cal- 
culated solution can in general not be accepted as fea- 
sible for the analytically posed problem — this means 
bad condition. As robust optimization yields a solu- 
tion which is feasible for a box including the original 
parameter and also the backward error caused by the 
backward error analysis performed for the stability of 
a specific algorithm, it provides a possible resource 
with artificially chosen limits. 



3 Robust Optimization 

The basic idea of robust optimization, as it will be 
used later, is outlined in the following. For the general 
optimization problem min^ z(x,y) subject to x E S(y) 
restricted to the fixed parameter y = yo with the ob- 
jective function z, which is without loss of general- 
ity real-valued, and the parameter-dependent search 
space S(y), a robust solution in the case of perturba- 
tions in the parameter y is searched. The perturba- 
tions constitute the set Y. Obviously, one has to dis- 
tinguish between objective function and constraints. 
In the sense of the worst-case-consideration only ro- 
bust feasible solutions X e S(y) := flyer S(y) are con- 
sidered. Clearly in general S{Y) could be empty and 
the problem is robust unfeasible. A robust optimal solu- 
tion is a robust feasible x which minimizes the worst- 
case objective function. With that follows the robust 
optimal solution subject to robust feasibility as solu- 
tion of minx^t i subject to Vy G Y : z(x,y) < f and 
X E S(y), cf. Ben-Tal and Nemirovski [1998]. In the 
linear case this robust optimization agrees also with 
Soyster [x ]. If y is not present in the constraints, 
then the robust optimization problem is equivalent to 
the well-known min-max problem originating from 
the game theory [Neumann and Morgenstern, 1944, 
chapter III]. 

4 Topology Optimization of Truss 
with regard to Volume (single 
load case) 

This general concept of robust optimization can now 
be applied to the topology optimization. 

Here, we do not treat the solvability. Instead we re- 
quest a reasonable problem and leave the question of 
the solvability to the user who provided it. 

The analytically stated and known (cf. [Przemie- 
iiiprVi^ 1 ^ matrix force method]; Marti [2003]) linear 
problem of the topology optimization of a truss with 
respect to its volume (plastic design) is given by the 
objective function to be minimized 

min Z^s (3) 

s,w 

with the bar lengths / E IR^ subject to the equilibrium 
condition 

Cw = f (4) 

with the reduced geometry matrix C G ]R^><^, the in- 
ner bar forces G and the reduced applied exter- 
nal loads / G IR^ . Further given are the constraints 

as <w< (5) 
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for the linear elasticity by — for all bars equal — stress 
limits > (7- G IR and < (7+ E IR (e.g. yield points 
for pressure and tension of the material all bars are 
made of), and the box constraint 

< S < Smax (6) 

for the cross sections of the bars s E K'^ with the given 
maximal bar cross section Smax E IR^. Here, a max- 
imum bar cross section adjusted to the bar length al- 
lows with (5) a prevention of buckling [Przemieniecki, 
1968]. The practical problem also includes the sup- 
ports able to absorb an arbitrary load. These support 
positions are removed from the equilibrium condition 
by canceling the respective lines, and it remains the 
reduced form (4) with the reduced loads and the re- 
duced geometry matrix. Here, we restrict ourselves 
to this reduced forms, and consequently abandon the 
explicit representation of the supports. 

In calculations relevant to the praxis, the stress lim- 
its (7- and (7+ are often reduced by a safety factor (i.e. 
0.7). The geometry matrix C describes the ground 
structure chosen in the modeling procedure. Hence, 
the only real parameter are the external loads /. As 
already shown in chapter 2, the topology calculated 
from the optimization problem (3)-(6) is bad condi- 
tioned with regard to the force direction. Therefore, 
now a substitute problem for (3)-(6) is sought, which 
is robust against perturbations / E f C IR^. 

For a robust feasible solution the constraints have to 



be fulfilled for all/G f: 

\ffeF:Cw{f)=f (7) 

(7_S < W{f) < (7+S (8) 

< S < Smax (9) 



Typically, the number of elements of F is infinite and 
therefore the number of constraints as well. Since the 
objective function (3) is independent of / and w{f) E 
R^, the robust optimization problem (3) subject to (7)- 
(9) arises. 

For a discrete approximation F = 
{firf2r---Jnp} C ]R^ of f with an at this point 
unspecified choice of follows the structure 

minfs (10) 

s,w 

s.t. V// E F : 

Cwi = fi E (11) 
as < Wi < (7+s (12) 

< S < Smax (13) 

of a linear program with the unknowns s E IR^ and 
Wi E K'^, also known as multiple load case. 

In the step to (10)-(13) we approximated the set F by 
F without stating how this step is performed techni- 
cally. This is in no case trivial! But manually choosing 



adequate discrete data from F to end up at F might be 
legitimate in a preliminary design, which is what the 
topology optimization often only provides. This man- 
ual choice is likewise not trivial. The perturbations 
yielding F have to be considered as a small neighbor- 
hood of /. A reasonable approximation F lies also in 
a small neighborhood of the applied /. However, by 
this naive approximation we were able to arrive at a 
computationally tractable problem (10)-(13) with the 
structure of the multiple load case. 

Nevertheless, with the technique presented in this 
paper one does not have to solve the problem of the 
approximation of the set F. In the following we are 
even able to treat a generalization of the single load 
case (3)-(6). 

5 Topology Optimization of Truss 
with regard to Volume (multiple 
load case) 

We now want to consider the multiple load case, 
which is apparently a natural extension of (3)-(6). In- 
stead of just one single load /, we look at a given but 

arbitrary set of loads f^/^ = |/i,/2/ • • • //n^ | C IR^. 
This leads to the optimization problem: 

min/^s (14) 

s.t. Mfi E f^/c : 

Cwi = fi (15) 
as < Wi < (7+s (16) 

< S < Smax (17) 

In contrast to f , the distance between two elements of 
F^lc cannot be assumed as small. We can again gener- 
ate a robust optimization against given perturbations 
F. With F^ic C F C follows as a substitution prob- 
lem with a typically infinite number of elements: 

min Z^s (18) 

s 

s.t. V/ E f : 

Cw{f)=f (19) 
as < w{f) < (7+s (20) 

< S < Smax (21) 

Assuming that f is a union of sets, i.e. F := 

[JiLi ^i' ^ill ^C)t need to approximate F, but rather 
we are able to specify an optimization problem which 
is equivalent to (18)-(21). Of course this procedure 
covers the special case of the single load case. The 
usual design of the limits as they are provided by the 
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engineer leads in the following to the assumption of F 
as the union of parallelotopes. The following is easily 
reproducible also for F as the union of general sets — 
but both the treatment and the solvability in a follow- 
ing optimization algorithm depend on the structure of 
these sets. 

Now, we examine the multiple load case which is 
given by the set f^/^ = |/i,/2, . . . | of ele- 
ments. Here, the force // acts in the nodal points given 
by the set Jj. We will describe the corresponding per- 
turbed forces // E Ff C by box constraints con- 
cerning the nodal points at which fi is acting: 



\fi-fi^fi + LE^(J'^)^i,iy' 

^ ie}ik=i 



A-(/,/,/c) < Si^j^k < A+(/,/,/c) 



(22) 



Here, d E {2,3} specifies the dimension of the 
space in which the structure lies. r(;,fc) G repre- 
sents a vector with the ;th component equal 1 for the 
;th node in the spatial direction fc, and elsewhere. 
A_(z,;,fc) E ]R and A+(z,;,fc) E K, respectively, repre- 
sent the perturbation in the ;th node for the zth force in 
the spatial direction k. In contrast to Ben-Tal and Ne- 
mirovski [1997], we consider the perturbations only 
in the nodes where the forces are acting. By choosing 
appropriate sets perturbations could also be con- 
sidered in any node, but here we want to obtain a ro- 
bustness with regard to the parameter force. 

To simplify the description, fj-^\ 1 < j < n/, 1 < / < 
fif are introduced as the edges of the parallelotope 
Fj for the z-th load case, so that the convex polytopes 
Fj can be written as the convex hull 

p. = convj/;^^^ : 1 <; < n,} 

of their edges. This allows to write every / E f/ as the 
convex combination 

7=1 7=1 

of the edges f-^'' for some appropriately chosen ocj. 
For every load case, we define S/ := 

|//^^ : 1 < ; < n/| as the set of the edges. Fur- 
thermore, we have 



Uf 



i=l z=l 

the set of all loads. 

To formulate the robust semi-infinite optimization 
problem (18)-(21) with a finite number of constraints, 
we need the following theorem: 



Theorem 1 Let 

Ml := {s E : V/ E f :3k;(/) : 

Cw{f) -- 
as < w 



/ 



and 



M2 := |s e R" : V/^- G (J S; Bwj 

i=l 



as 



< Wj < (7+5 1 . 



(23) Then 



Ml = M2. 



Proof We have to show Mi C M2 and M2 C Mi. 

In the first case Mi C M2, let s E Mi be arbitrary. 
For an arbitrary /y E U/£i holds: 

fj E F ^1^^: 3w{fj) : as < w{fj) < a^s 

Since /y was arbitrary, for every /y exists ^{fj). This 
means s E M2. 

To proof the second case M2 ^ Mi, we choose an 
arbitrary s E M2. Let / E f , then there exists 1 < z < 
n f with f e Fj = conv Sj , and we can write / as the 
convex combination 



(7) 



f = L-ifi 

7=1 



for appropriate a:i, a:2, • • • , a:^^ with Y^^^-^ocj = 1. As 
s E M2, there exists Wj: 



3wi,...,Wn^ :Cwj = f^^"^ 

as < Wj < a^s 

So that for w{f) := Ejii ocj'^j holds 



Cw{f) = C ^ A:yK;y 
7=1 



(24) 
(25) 



(26) 
(27) 
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and 



W = OCjWj 
7=1 

(25) ( \ 

y=i Vy=i / 

W = ^ OCjWj 
7=1 
(25) ^ 



(28) 



(29) 



(30) 



(31) 



Since / was arbitrary, (26)-(31) holds for all / G P. 
This means s G Mi . ■ 

With S := |J/!£i 5/ we can now reformulate the semi- 
infinite optimization problem (18)-(21) with a finite 
number of constraints: 



min/^s 

s,w 



s.t. \/fj e S : 

Cwj = fj 

as < Wj < a+s 

^ S < Sjnax 



(32) 

(33) 
(34) 
(35) 



To get a typical linear program, we consider the im- 
plicit variables w as design variables. Theorem 1 en- 
sures equal polyhedrons in (18)-(21) and (32)-(35) and 
thus theses linear programs are equivalent. 

Considering a more general F as the union of con- 
vex sets Fj, the convex hull of the extreme points — in 
a linear case the extreme points are the edges — de- 
scribes the same sets. Therefore, for finite numbers of 
extreme points of the sets the above approach can be 
performed in the same way. 



is pS with the same topology. The equation Cw = af 

is a short hand for Cw = / with a scaled design space 
represented by the scaled geometry matrix C = ociC, 
a scaled force / = 0.2/ and a: = ^, < oc,oci,oc2 G K. 
The same holds for the inequation. 

The cross section size of all bars is scaled for the vi- 
sualization, but in Figs. 1, 2, 3, 4, 5 and 6 the same 
scaling is used. The Figs. 7 and 8 have the same scal- 
ing as to each other, too. 

The calculation is done by the library from 
Makhorin [2010]. We used the simplex solver as well 
as the interior point solver. 

Example 1: In the following first example (c.f. Figs. 
1, 2), the maximum quantity of all loads is equal 10^ N. 
For some chosen values, the initial examples result 
in the numerical solutions s^^^^ for a strictly horizon- 
tal force (Fig. 1(b)) and s^'^\ shown in Fig. 2(a), for 
fy = 0.1/.: 



0.00 
1.00 • 10" 
0.00 



s(2) ^ 



1.41 • 10-5 
8.96 • 10-5 
0.00 



For a symmetric perturbation of the horizontal force 
of 10 % in every direction the solution of (32)-(35) is 
the robust solution s^'^\ illustrated in Fig. 2(b): 

/ 7.81 • 10-°*' 
1.21 10-°^ 
\ 7.81 • 10-°"^ 



c(3) 



The optimal values for the objective function are the 
volumes in m^: 



optimal 


perturbed 


robust 


0.000100 


0.000109 


0.000144 



6 Numerical Examples 

In all examples we use as material an aluminium with 

a yield strength of 10^ Pa and a density of 2.7 • 10^ 
The measurement unit for the sizes of the design 
spaces is meter. 

All these data are only scaling factors in the linear 
program and do not affect the resulting topology. For 
example let us consider s as a solution of: 

min/^s s.t.: Cw = f;(J-S <w< c7+s;0 < s 

s,w 

Then for < a:, j6 G IR a solution of 

min/^s s.t.: Cw = ocf;BcT-S < w < Bct^s;0 < s 

s,w 



Example 2: The prime example of a truss from 
Ben-Tal and Nemirovski [1997], there calculated by 
semi-definite programming with regard to robustness 
based on a finite element approach, is illustrated in the 
figures 3(a)-3(c), but now calculated with our robust 
optimization method. Here, the maximal quantities 
of the perturbed forces are equal to the unperturbed 
forces. 

Figure 3(a) shows the ground structure with the 
fixed nodes on the left side. The arrows repre- 
sent the forces of equal quantities of 10^ N. Figure 
3(b) displays the optimal solution with a volume of 
0.000800 m^. The robust optimal solution with a vol- 
ume of 0.001309 m^ for a symmetric perturbation of 
10 % in every direction of every force is presented in 
figure 3(c). 
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Figure 3: (a) Ground structure of example 2 with fixed nodes on the left side and arrows representing the 
forces, (b) Optimal solution of the example 2 shown in Fig. 3(a). (c) Robust optimal solution of example 2 (c.f. 
Fig. 3(a)). 



The results of the topologies are congruent to the 
results in Ben-Tal and Nemirovski [1997]. Also, in 
Calafiore and Dabbene [2008] the same topology is ob- 
tained by a sampling-based approximate solution of 
the worst-case. 

Example 3: We calculated a 3-dimensional example 
with 27 nodes {a,b,c) e K^, a,b,c G {1,2,3} mea- 
sured in meter and 274 potential bars between every 
two not fixed nodes — long bars which are located 
alongside several shorter bars are ignored. The sup- 
ports are in the nodes {l,a,h), a, b E {1,2,3} which 
are fixed in every direction. In the node (3,2, 1) acts a 
force of 4 • lO'* N in the negative z-direction. The opti- 
mal solution with a volume of 0.0024 m^ is presented 
in Fig. 4. For a perturbation of 10 % in every direction 
of the force, and scaling to obtain the same maximal 
quantities as in the unperturbed case, our robust op- 
timization method produces the robust optimal solu- 
tion shown in Fig. 5 with a volume of 0.0026 m^. 

Both Fig. 4 and Fig. 5 show two solutions which are 
caused by linear programming. The calculation with a 
simplex algorithm results in edges of the polyhedron 
shown in Fig. 4(a) and 5(a), whereas the results shown 
in Fig. 4(b) and 5(b) are calculated by an interior point 
algorithm which results typically not in an edge of the 
polyhedron. Hence, every point on the segment be- 
tween the points represented by Fig. 4(a) and 4(b), 
and Fig. 5(a) and 5(b), respectively, is also a solution 
with the same objective function value. 

Example 4: In this example of a mast are 72 nodes 
(a,b,c) E ]R^a,b E {0,i,^},c G {f : x = 0,1,. ..,7} 
measured in meter, and 503 potential bars between 
every two not fixed nodes with a maximal distance 

of ^ — long bars which are located alongside sev- 
eral shorter bars are ignored. The supports are in the 
nodes (a, b,0), a, b G {0, ^}, which are fixed in every 

direction. In the node ^, 1^ acts a force of 4 • 10^ N 

in the negative z-direction. The optimal solution with 




Figure 4: Optimal solutions of example 3. (a) Simplex 
algorithm, (b) Interior point algorithm. 




Figure 5: Robust optimal solutions of example 3. (a) 
Simplex algorithm, (b) Interior point algorithm. 
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(a) 



(b) 



Figure 6: (a) Optimal solution of example 4. (b) Robust 
optimal solution of example 4. 




Figure 7: (a) Ground structure of example 4. (b) Opti- 
mal solution of example 4. (c) Robust optimal solution 
of example 3. 



a volume of 0.000514 m^ is presented in Fig. 6(a), 7(b). 
For a perturbation of 50% in every direction of the 
force and scaling to have the same maximal quantities 
as in the unperturbed case, our robust optimization 
method produces the robust optimal solution shown 
in Fig. 6(b), 7(c) with a volume of 0.001568 m^. 

Example 5: As a last example, we show a multiple 
load case. It represents the top part of a transmis- 
sion tower which supports two ground wires. The 

27 nodes (a,b,c) G K^, a,b,c e |o, ^/l} are con- 
nected by 294 potential bars between every two not 
fixed nodes — long bars which are located along- 
side several shorter bars are ignored. The 4 nodes 

(i,0,0),(i,l,0),(0,i,0),(l,i,0) G ]R3 are fixed 

in every direction. We have 2 forces fi and /2 act- 
ing alone or together constituting the 3 load cases — 
represents the existing of either each single or both 
ground wires on top of the transmission tower, fi and 
/2 are acting in (0, 0, 1) and (1, 1, 1), respectively, in the 
negative z-direction with a magnitude of 2 • 10^ N. 

The optimal solution with a volume of 0.000850 m^ 
is presented in Fig. 8(a). For a perturbation of 50 % in 



every direction of the forces, and scaling to obtain the 
same maximal quantities as in the unperturbed case, 
our robust optimization method produces the robust 
optimal solution shown in Fig. 8(b), 8(c) with a vol- 
ume of 0.002562 m^. Fig. 8(b) is the result of a simplex 
algorithm and therefore an edge of the feasible region, 
whereas Fig. 8(c) is calculated by an interior point al- 
gorithm. 

7 Conclusion 

In the introduction we became acquainted with the 
origins of uncertainties. In our opinion every treat- 
ment of uncertainties leads in the end to robust opti- 
mization. The bad condition of the topology of a truss 
makes it necessary to choose artificial limits to defi- 
nitely get a feasible solution for the analytical prob- 
m from the computer. In contrast to the common fi- 
ite element approach of topology optimization, the 
asic initial approach chosen here is able to mini- 
mize the volume. The robust optimization yields a 
semi-infinite optimization problem. Typical commer- 
cial software products are able to handle semi-infinite 
optimization problems by adaptive generation of fi- 
nite approximations of the infinite constraints. The ex- 
pert knowledge of an engineer is sometimes also able 
to generate a wise approximation. Every time, these 
discrete approximations lead to the structure of the 
multiple load case. If we are looking for a black-box 
computer program to solve our problems, we need a 
well-defined and straightforward procedure. Our ap- 
proach is able to fulfill this dream without an approx- 
imation. Finally, the numerical examples show the 
advantage of the robust optimization with the only 
drawback of a slight increase of the volume. 
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